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Abstract. In this work we study the properties of a relativistic mixture of two non-reacting species 
in thermal local equilibrium. We use the full Boltzmann equation (BE) to find the general balance 
equations. Following conventional ideas in kinetic theory, we use the concept of chaotic velocity. 
This is a novel approach to the problem. The resulting equations will be the starting point of the 
calculation exhibiting the correct thermodynamic forces and the corresponding fluxes; these results 
will be published elsewhere. 
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INTRODUCTION 

The purpose of this work is to obtain the balance equations for a system composed of 
two inert, non degenerate, dilute species. The calculation will be carried out within the 
framework of special relativity. Its importance today relies on the fact that many appli- 
cations can be found in a variety of problems, from astrophysical systems to relativistic 
heavy ion collisions. Contrary to most of the approaches which are given in the literature 
for relativistic gases we here introduce, as in classical kinetic theory, the concept of 
chaotic velocity. This will be reflected on the fact that the resulting fluxes will be 
analogue to those obtained more than one hundred years ago by Maxwell and Clausius 
[1, 2, 3], for whom the concept of chaotic fluxes was key in describing dissipative effects. 

In particular, we here address the first stage of this program namely, the derivation 
of the balance equations for the state variables. In our case these are chosen to be, 
the number densities of both species, the baricentric velocity of the mixture and the 
internal energy density. Further the fluxes will be calculated in an arbitrary frame as 
the transformed dissipative fluxes which are identified in the local co-moving frame. 
This choice will reflect the importance of the Lorentz transformations in identifying the 
concept of chaotic velocity. 



TRANSPORT EQUATIONS FOR THE NON-RELATIVISTIC 

BINARY MIXTURE 



Kinetic theory constitutes a microscopic formalism from which the macroscopic prop- 
erties of a system can be obtained based on the evolution of a distribution function for 
the molecules in the gas. For a binary mixture one considers two separate distribution 
functions, one for each species, which satisfy two coupled Boltzmann equations, namely 

^ + Va~=J(fafa)+J(faf b ) (D 

where the collision term is given by 

J(fifj) = /-/ [fWfiv'j) ~ fWfivj)] x a {y{vj yWj) gij^MK (2) 

for each species ij = a,b. Notice that the coupling of both integrodifferential equa- 
tions is given by the second term on the right hand side of Eq. (1) namely, the cross- 
collision term. The primes denote the values of v, after the binary collision takes 

place, o (viVj — > v | v yj dvjdv^- is the cross section, namely, the number of molecules 

per unit time of species i colliding with a molecule of species j such that after the 
collision the molecules have velocities v- in the range dv- and v^- in the range dv-; 
gij = |vj — Vj| = |v- — y'j\. We also recall the reader that the cross section o satisfies 
the principle of microscopic reversibility, namely, it is invariant upon spatial and tempo- 
ral reflections, so that, 

a (v,-v j ->• v-y'j) = o (v-v^- ->■ v ; -v j) for i,j = a,b (3) 

thus guaranteeing the existence of inverse collisions. 

Since the treatment in both kinetic equations is symmetric we will only refer to species 
a. There is an analogous procedure for species b. In order to obtain the balance equations, 
the solution of the Boltzmann equation is not required since it already contains the 
information about the conservation of macroscopic quantities through the invariance of 
mass, momentum and energy in each individual collision. As usual, we define the local 
particle density as, 



n a = J /adv a , 



(4) 



To obtain the conservation equations, one multiplies Eq. (1) by the collision invariants 
and integrates over velocity space. For the mass conservation, we multiply by 1 and 
obtain 

^ + V.(p fl u) = -V.J fl (5) 



where p a = m a n a such that p = p a + p b is the total density. The baricentric velocity u is 
defined as 



mi = £n/Uj = / v */i dv n ( 6 ) 

i'=a i=a 

and the mass flux, 

ia = m a J k a f a dk a , (7) 

where the chaotic velocity for species a is given by 

k a = \ a - u, (8) 

and the flux satisfies 

ia = ~h- (9) 



For the momentum balance equation, we consider the collisional invariant m a \ a . With 
the help of k a = y a — u, one obtains 

^(pu) = -V-(T + puu), (10) 



where 



b f 

T=Y d m i jk i k i f i dk i (11) 



is the stress tensor. 



Finally for the energy balance equation we use the invariance of \m a \y a \ 2 , and in a 
similar fashion we obtain, 

— (pe) = -V-(peu + J ? )-T:Vu, (12) 

where 

b b m, 

Hi 



i=a i=a 



J<7 = £j<7,- = £- Wifidki, (13) 



is the heat flux, and the local internal energy is given by 

b 



pe = ^m l jkff l dk l , (14) 



where k[ is the magnitude of k,-. 

It is important to emphasize the fact that the averages in Eqs. (7), (11) and (13) are 
calculated with the chaotic part of the molecular velocity k a . Isolating the chaotic part of 
\ a leads to the identification of the dissipative fluxes. In fact, his is the physical meaning 
of heat firstly introduced by Maxwell [1] [3]. 



RELATIVISTIC KINETIC THEORY IN SPECIAL RELATIVITY 

The system we consider, and that will refer to as relativistic binary mixture, is a gas 
constituted of two non-degenerate species that do not react but only interact through 
collisions. The gas is diluted but the molecular velocity of both species is high enough 
for relativistic effects to be relevant. This is reflected in the fact that the relativistic 
parameters zt = kT/mic 2 are grater than one namely, the thermal energy is larger than 
the rest energy of the molecules. Here T is the temperature of the gas, k the Boltzmann 
constant, c the speed of light and m,- the rest mass of each species. Additionally, we 
consider the system in the absence of external forces. Before addressing the kinetic 
theory of the mixture we review some basic aspects of the relativistic kinetic theory. The 
relativistic Boltzmann equation reads [16, 17, 18], 

v a f,a=J(ff) (15) 

where / is the distribution function as before. Proofs of the invariance of this quantity 
are available in Refs. [17] and [18]. Here v a = y ffl ( G^,c) is the molecular four-velocity 

I Irs 

with 7o= (l — (co/c) 2 ) being the Lorentz factor. Thus, the left hand side of Eq. 
(15) is clearly an invariant. The collision term can also be written in an invariant fashion 
as follows 

J(ff) = J J {ffi-ffi)*o(a)dadv* u (16) 

where dvT = ^-p- and & are Lorentz invariants. The latter is known as the invariant flux 
and is given by [18], 

& = ^V\\ = ~ c \/^V lv ) 2 -C 2 = ^VW»1 («1- C 2)) 2 - C 4 (17) 

which reduces to the relative velocity in the non-relativistic limit. 

It is well known that without a solution of Boltzmann equation we can obtain two 
results, the H theorem and the balance equations. Since in this work we focus on the 
balance equations we will not discuss any further the properties of the Boltzmann 
equation or its methods of solution. 



As in the non-relativistic case, the Boltzmann equation is multiplied by each of 
the collision invariants and integrated to yield the transport equations. However, the 



dissipative fluxes cannot be clearly identified if the chaotic velocity is not introduced 
explicitly. As extensively discussed in Ref. [21], Lorentz transformations can be used in 
order to write the molecular velocity measured by an arbitrary observer in terms of the 
chaotic velocity which is the one measured in a local frame that moves with the fluid 
element namely, the co-moving frame. This is expressed as 

v^=^K v (18) 

where K v = fk (jt , cj is the chaotic velocity. The establishment of the particle four-flux 

and energy-momentum tensor in this context for the one-component system is discussed 
in Ref. [21]. In the next section we will follow these ideas in order to obtain the balance 
equations for the relativistic mixture. 

BALANCE EQUATIONS FOR THE RELATIVISTIC MIXTURE 

The Boltzmann equations for the special relativistic binary mixture are 

vf/i,« = /(/i/i)+/(/i/ 2 ) (19) 

V?/2,«=/(/2/2)+/(/2/l) (20) 

where J(fifj) is defined in the same way as in Eq (16) and the indices 1 and 2 indicate 
species. 

As before, in order to obtain the balance equations one multiplies Eqs. (19) and (20) 
by the corresponding collision invariants, in this case m, and vf, i = 1,2 and integrates 
over velocity space. By multiplying Eq. (19) by m\ and integrating over the Lorentz 
invariant element dv\ one finds, 

<u = (21) 

where 

N» = < +N$ = J vf/idvj + J v»f 2 dvl (22) 

which corresponds to the four-flux of particles in an arbitrary frame. Equation (22) 
leads to the definition of a baricentric velocity given by W =nll^, where n = n\ + n 2 . 
Notice the similarity with the definition in the non relativistic case, Eq. (6). By using the 
transformation in Eq. (18) one can write 

N? = 2fjf, (23) 

where J* is the dissipative particle four-flux in the co-moving frame and it satisfies 

J\ -(-/J = (0,^)- Introducing this relation in Eq. (21) one finds after some laborious 
algebra, that 



dt 



J i + ^ uUi §-/i +n l d + U^ = 0, (24) 



c 



where U v = y u («,c), y u = (l — (u/c) 2 ) and = t/^. Equation (24), which corre- 
sponds to the mass conservation equation for species 1; latin indices correspond to the 
spatial part of any tensor and run form 1 to 3. Notice that in Euler's regime only the last 
two terms do not vanish. 

For the energy momentum conservation we multiply Eq. (19) by vj\ integrate on 
velocity space and obtain, 



7?;=o 



where 



T» v = Tr + T? v = mi J v f tv\f 1 dvt + m 2 J vf^v v 2 f 2 dv* 2 , 
is the energy-momentum tensor. 



(25) 



(26) 



In order to explicitly calculate Eq. (26) we start by noticing that the energy- 
momentum tensor, in the co-moving frame can be written as, 



( p \ 

p 

p 

V e / 



+ 



^11 ^12 ^13 <?1 \ 

^12 ^22 ^23 12 

^13 K 23 K 33 <?3 

q\ q 2 q3 K44 J 



(27) 



where first term corresponds to the equilibrium situation [20]. Here the hydrostatic 
pressure p and internal energy density £ are given by 



P = Pl+P2 = m l J KfK^d 3 ^ +m 2 j K^ff ] d 3 K : 



(28) 



e = niei + n 2 e 2 = mi J « f^d 3 ^ + m 2 J K^f^d 3 ^, (29) 

where denotes the local equilibrium solution to Eq. (19). This function is the well 
known Jiittner distribution function [19],[17][18]. 

The second term in Eq. (27) corresponds to the non equilibrium situation, where we 
identify, 



q a = q a + q a = c 2 J y^f^R* + c * J y^frdK^ 



(30) 



and 



7t ab = Jt ab + n ab 



m x J K?K b 1 f l dK* 1 +m 2 J K?K%f 2 dK}. (31) 



In Ref. [21] it is shown that 



n 44 = (32) 

by introducing a Chapman and Enskog expansion. The quantities in Eq. (30) and (31) 
are identified as the heat flux and viscous tensor since they correspond to the average 
of the chaotic energy and momentum fluxes respectively. This is based on the physical 
interpretation of dissipative fluxes as is mentioned in Refs. [1][2][3]. 

We then proceed in the same fashion as in the particle conservation, by writing 

T^ v = &£&£T a t (33) 

which yields 

T a $ = pg a P + 1 (p + e) U a UP + 1 (u a ^q a + U^ a a q a ^j + ^^ b a Tl ab , (34) 

the first two terms correspond to the equilibrium case, while the second term represents 
the contribution of the heat flux and the third term includes the viscosities. Here 
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The last step is to work out the derivative in Eq. (25) with the help of Eq. (34). This 
leads to 



1 



U P eU^n^+nU^U^e^+neU P U^ + neU^U^ 



+P, a h^ + ^(UPU V V + U^) (35) 
U a J?f) f + (Wj 3 ) q% + (V^f ) J + (UPJK?) q% 

n ab + (j2?js?* a ) K = o, 



which is the energy momentum balance equation for the mixture. In Euler's regime it 
reduces to, 



[ne + pO] + \ (ne + p)U l5 + p. a h^ a = 0, 



(36) 



where the notation ( ) = U^( ) ;jU and = has been introduced. To isolate the energy 
balance, we calculate UpT.^ v which, again in Euler's regime, reads 



(ne + p6) = 0. 



(37) 



Introducing the previous equation in Eq. (36) we obtain 

\(ne + p)U p +h pv p }V = 0. (38) 

It is important stress that the form of Eq. (34) has been obtained form microscopic 
grounds. 

CONCLUSIONS 

We have obtained the balance equations for a relativistic, diluted, non degenerate, 
mixture based solely on microscopic grounds. We explicitly used the concept of chaotic 
velocity to identify the dissipative fluxes. The results are shown in Eqs. (24) and (35) 
from which the Euler equations are recovered in the equilibrium case. The dissipative 
terms in these equations differ from the ones obtained following the standard procedure 
(see appendix B of Ref. [21]). 

In particular, the Euler equations for the mixture are required in order to solve 
the linearized Boltzmann equation within the Chapman and Enskog expansion. This 
constitutes work in progress and will be publish elsewhere. 
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